Correlated time-dependent transport through a 2D quantum structure 
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We use a generalized master equation (GME) to describe the nonequilibrium magnetotransport 
of interacting electrons through a broad finite quantum wire with an embedded ring structure. The 
finite quantum wire is weakly coupled to two broad leads acting as reservoirs of electrons. The mutual 
Coulomb interaction of the electrons is described using a configuration interaction method for the 
many-electron states of the central system. We report some non-trivial interaction effects both at 
the level of time-dependent filling of states and on the time-dependent transport. We find that the 
Coulomb interaction in this non-trivial geometry can enhance the correlation of electronic states in 
the system and facilitate it's charging in certain circumstances in the weak coupling limit appropriate 
for the GME. In addition, we find oscillations in the current in the leads due to the correlations 
oscillations caused by the switched-on lead- system coupling. The oscillations are influenced and 
can be enhanced by the external magnetic field and the Coulomb interaction. 

PACS numbers: 73.23.Hk, 85.35.Ds, 85.35.Be, 73.21. La 



I. INTRODUCTION 

The Coulomb blocking of electrons entering nanostruc- 
tures has been known for quite some time and is ex- 
plained by the magnitude of the direct part of the re- 
pulsive Coulomb interaction energy in relation to the 
energy spectrum of the nanosystem.^ The experimental 
and theoretical studies of Coulomb effects on the meso- 
scopic transport were commonly focused on steady-state 
regime of single or double quantum dots. However, the 
increasing interest in fast dynamics at nanoscale and 
time-resolved detection of electrons via a nearby detector 
strongly motivates theoretical investigations of interact- 
ing time-dependent transport in complex systems. The 
dynamical aspects of the Coulomb blocking have been 
investigated by Kurth et al. in a one-dimensional lat- 
tice model using combination of non-equilibrium Greens 
functions and time dependent density-functional theory 
for the Coulomb interactionP 

In a recent work 3 we have also analyzed the transient 
currents through an interacting two-dimensional quan- 
tum dot by solving the generalized master equation for 
the matrix elements of the reduced density operator act- 
ing in the Fock space of interacting many-electrons states 
of the dot. The GME scheme that we implement numer- 
ically takes into account the geometrical details of the 
sample and leads (see Ref. [4]). Our aim here is to ap- 
ply the same method to more complex systems in order 
to single out non-trivial Coulomb effects in the transient 
regime. The system we consider is a parabolic quantum 
wire with an additional ring-shaped confining potential. 
The embedded ring geometry imposes different localiza- 
tion properties of the states with respect to the regions 
where the leads are attached. This fact has important 
consequences on the time-dependent filling of the many- 
electron states. 



Even though the transport and magnetic properties of 
quantum rings have fascinated researchers for a long time 
new questions and results regarding the Aharonov-Bohm 
interference oscillations have been catching attentionP^ 
Dynamical effects have been studied in open or closed 
quantum rings. Propagation of electron pulses in rings 
of finite width has been investigated by Chaves et alPand 
by Thorgillsson et alP within scattering theory, and non- 
adiabatic current generation in a closed finite quantum 
ring in an external magnetic field has been studied by 
integrating the Liouville-von Neumann equation for the 
density operator in timeP^The system was perturbed by 
a strong dipole or higher order multipolar electric electric 
field pulse and the mutual Coulomb interaction between 
the electrons was included in a mean-field DFT manner. 

We use a non-Markovian version of the generalized 
master equatiorPEl anc [ treat the Coulomb interaction of 
the electrons exactly within a truncated many-electron 
basisP The GME formali sm w as originally proposed 
by Nakajima and Zwanzigj 12 * 13 * and has more recently 
been applied to study transport phenomena by several 
authors j 14 " 17 ! just to cite few. 

The paper is organized as follows: Section II sets the 
notations and presents briefly the formalism, Section III 
describes the quantum wire structure with the embedded 
ring, Section IV contains the results and their discussion 
while Section V is left for conclusions. 



II. THE GME AND THE COULOMB 
INTERACTION 

In this Section we recall the main outlines of the GME 
method leading to the numerical results. The time- 
dependent transport of noninteracting electrons with 
the generalized master equation has been described in 
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two publication for a lattice model 4 and a continuous 
model. 11 The non-interacting many-electron Hamilto- 
nian of the coupled system (i.e. the central sample, the 
leads and time-dependent coupling) reads as 

H(t) = J2E a dUa+ £ e l (q)cl lCql + H T (t), (1) 

a q,l=L,R 

where the tunneling Hamiltonian i^x(^) = H^(t)-\-H^ (t) 
describes the coupling of the system to the left and right 
leads 

H l T (t) = X \t) {T l q J ql d a + {T l qa )*dlc ql ). (2) 

q,a 

The creation and annihilation operators are associated 
to single particle states of the disconnected subsystems. 
The coupling coefficient T l of a single-electron state \q) 
in the lead / to a state \a) in the system is modelled as 
a non-local overlap integral of the corresponding wave 
functions in the contact regions of the system, Q l s , and 
the lead I, 

T l aq = [ dr^(^(rO)*^(r)4(r,rO + ^c. (3) 

Jn l s xni 

The function 

«4(r,r') = sfcexp [-5[(x - x'f - 5 l 2 (y - y'f] 

with r G and r' G £li defines the 'nonlocal overlap' and 
their affinity in energy. The semi-infinite leads have the 
same parabolic confinement as the finite quantum wire 
in the ^/-direction, perpendicular to the transport direc- 
tion x. The confinement is characterized by the energy 
scale hQo, The energy spectrum of the leads e l (q) is con- 
tinuous, but with clear subband structure. The effects 
of the external magnetic field B = Bz is present in the 
energy spectrum of the leads e l (q), the spectrum of the 
system E a , and in the wave functions of the leads and 
the system. 

In order to describe the time-dependent transport 
when the system contains few electrons we select the 
lowest A/ses single electron states (SES) of the central 
system to construct a Fock space with TVmes = 2 NsES 
many-electron states (MES). In the occupation represen- 
tation basis such a state can be written as 

| /i } = K,^,...,C---), (5) 

where is the occupation of the n — th single parti- 
cle state of the isolated system. A/ses is selected large 
enough that the chemical potentials of the leads \i\ in 
equilibrium before the coupling at t = are smaller than 
the energy of the highest SES, and ideally a further in- 
crease of A^ses should not change the transport results 
of the calculations. 



The Liouville-von Neumann equation describing the 
time-evolution of the total system, central system and 
leads 

ihW(t) = [H(t), W(t)], W(t < t ) = PlPrPs, (6) 

where the equilibrium density operator of the discon- 
nected lead I with chemical potential \i\ is 

^(Hi-fuNi) 

pl = Tr z {e-^-^^)}' 

is now projected on the central system by partial trac- 
ing operations with respect to the operators of the leads. 
Defining the reduced density operator (RDO) of the cen- 
tral system 

p(t) = Tr L Tr R W(t), p(t ) = ps, (8) 

we obtain an integro- differential equation for the RDO, 
the generalized master equation (GME) 

Pit) = ~[H s ,p(t)] 

- i E / dq X l (t)([T l Mt)] + h.c), (9) 

where two operators have been introduced to compactify 
the notation 

J t 

IW = Us(s) (T'V(s)(l - /') - P {s)T^f) Ul(t), 

with Us(t) = e*^/ a ^ s , and a scattering or coupling op- 
erator T acting in the Fock space of the system 

T\q) = Y, T UlM)((3\ (10) 
TUq) = £lt(a|4|/3). (11) 

a 

Here the kernel of the integro-differential equation has 
been obtained by taking into account only second order 
processes with respect to the coupling coefficients. It 
should though be kept in mind that the structure of the 
equation implies higher order processes to infinite order. 

In the derivation of the GME here only the coupling 
Hamiltonian is allowed to depend on time. The possibil- 
ity of the Hamiltonian of the central system Hs depend- 
ing on time (describing a laser pulse for example) has 
been considered by Amin et al. 18 A DFT description of 
the mutual Coulomb interaction of the electrons in the 
central system would require Hs to be time-dependent, 
and it would also require a further reduction of the GME 
introducing the reduced single-particle density matrix 
loosing some many-electron correlation effects caused by 
the coupling of the system to the leads 
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Our approach to solve this dilemma has been reported 
earlier pi but here we shall briefly outline it for the case of 
the continuous model. We choose to change the Hamil- 
tonian of the central system 

H s = Y,E a d\d a + l - Y,{ab\V C oux\cd)d{d\d d d c (12) 

a abed 

to include the time-independent Coulomb interaction 
term appropriate for our many-electron formalism. The 
Coulomb potential is 



^Coui(r - r') 



KyJ (x — x') 2 + (y — y') 2 + r\ 2 
and the matrix elements are 



, (13) 



(ab\V Co ui\cd) = / dr*:(r)/ 6c (r)* d (r) (14) 



with 



4c(r) 



J dr'^* b (r')Vcou^r-r')^ c (r'), 



(15) 



and r] is a small convergence parameter to be specified 
later. Along the lines of approaches developed under the 
names of "configuration interaction" or "exact numeri- 
cal diagonalization" we diagonalize the new interacting 
Hamiltonian ( [l2| ) in the MES basis of the non-interacting 
system {\p}} in the entire Fock space built from the A$es 
SES states, including all sectors containing zero electrons 
(the vacuum state) to Ases electrons, since we are deal- 
ing with an open system with variable number of elec- 
trons. The diagonalization yields a new basis of interact- 
ing MES {\p)} connected to the non-interacting one by 
a unitary transformation 



Im) = 5Z v ^l a )' 



(16) 



supplied by the diagonalization. Here we need to keep 
in mind that V will be represented by an Ames x ^mes 
matrix in numerical calculations. An inspection of the 
structure of the non-interacting GME ([9| reveals that the 
equation can be also be transformed to the interacting 
basis {|/i)} by the unitary transformation. Thus, in a 
numerical calculation a basis transformation of the many- 
electron coupling matrix T l (q) = V^T l (q)V (11) and the 
insertion of the diagonalized matrix representation of the 
interacting Hs in the GME (|9| will give the RDO in 
the interacting MES basis p = V^pV. As all measurable 
quantities are in the end expressed as a partial trace with 
respect to operators of the central system expectation 
values can be calculated in the new basis and for the same 
reason in the non-interacting case we can again obtain the 
mean value of the left or right current directly from the 
transformed GME. 

In our earlier publication^ neglecting the Coulomb 
interaction it was clear that we could only effectively 



describe the time-dependent transport through systems 
with up to 7 or 8 SES considered relevant for the cur- 
rents, situated in and around the window of the chemical 
potentials of the two leads. This limitation was imposed 
by the complex bandstructure of the energy spectrum 
in the broad leads employed in the calculations. Here, 
we may have to include more SESs in order to describe 
reasonably the interaction of the electrons in the central 
system. In order to accomplish this we have to resort 
to a more refined truncation procedure than we used for 
the non-interacting system: The unitary transformation 
can not be truncated and has to include the Ames states 
constructed initially from the Ases SESs. In the numer- 
ical calculations here we will employ 12 SES leading to 
4096 MES. The unitary transformation of T l (q) is thus 
CPU-time intensive for all the q- values necessary for the 
leads, but it has only to be performed once. After that it 
is possible to deploy a second truncation to the GME by 
keeping only the A" MES << Ames MES with lowest en- 
ergy. Typically, for the parameters that we will select for 
the numerical calculations here we need only A^ MES = 32, 
but this cutoff is very system dependent. 



III. EMBEDDED QUANTUM STRUCTURE 

Here we will use the GME to analyze time-dependent 
transport of electrons through a short but broad quantum 
wire of length L x = 300 nm, with an embedded quantum 
ring. The parabolic confinement of the quantum wire is 
characterized by the energy scale hQo = 1.0 meV, and we 
assume GaAs parameters with m* = 0.067m e and n = 
12.4. The embedded ring is represented by the potential 



^QR(r) 



; Vi exp [-{p xl xf - (f3 yt y) 2 } + V g , (17) 



i=l 



with Vi = -4.0 meV, V 2 = +14.0 meV, f3 xl = 1.09 xl0~ 2 
nm" 1 , /3 y i = 3.46 x 10" 4 nm -1 , f3 x2 = 1.09 x 10" 2 nm" 1 , 
and f3 y 2 = 2.83 x 10 -2 nm -1 . The parameter V g can be 
thought of as a gate voltage. We use it to position the 
chemical potential of the right lead pr at a similar place 
in the energy spectrum of the SESs for the two values 
of the magnetic field investigated here. (For B = 1.0 T 
we use V g = 1.0 meV and for B = 0.5 T we have V g 
1 .2 meV). Figure ([!]) shows the ring embedded in the 
quantum wire with the spatial coordinates scaled by the 
magnetic length modified by the parabolic confinement 
a w = y/h/(m*n w ), with £l 2 w = ^ + uj 2 c at B = 1.0 T, 
where the cyclotron frequency is uo c = eB/(m*c). At 
B = 1.0 T a w = 23.87 nm. In the following calculations 
we use for the coupling (H) S[ = 4.39 x 10 -4 nm -2 , = 

/ 3 /2 

0.5 meV, and g dw = 30 meV or 40 meV. 

The energy spectra for the closed system of a quan- 
tum wire with an embedded ring are shown for the SESs 
and the MESs in Figure [2] with the chemical potential 
of the right lead to be used in the following dynamical 
calculations pr indicated. As will be evident from the 
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FIG. 1: (Color online) The potential denning the quantum 
ring (17) embedded in the finite quantum wire at B = 1.0 T 
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tively dense spectrum of MESs occupied by 2 electrons. 
Due to the strong Coulomb interaction the lowest MES 
occupied by 3 electrons is located well above the highest 
energy shown in Fig. [2] 

The value of hr = 1.6 meV is selected such that below 
it there are both SESs localized away from the contact 
region and states with a strong weight in that region. 
The probability for the SESs is displayed in Figure [3J 
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FIG. 3: (Color online) The probability density of the single- 
electron eigenstates of the system labeled by a in numerical 
order with a = 1 at the top left and a = 12 at the bottom 
right. B = 1.0 T, L x = 300 nm. 



IV. DYNAMICAL TRANSPORT PROPERTIES 



FIG. 2: (Color online) The energy spectra for the single elec- 
tron states (SESs) and the many-electron states (MESs) built 
from the lowest 12 SESs. The solid green line indicates the 
chemical potential of the right lead [1r = 1.6 meV. B = 1.0 
T, V g = 1.0 meV, L x 300 nm. 

probabilities of the SESs the two lowest SES states are 
almost degenerate, so below jiR there are 6 SESs. In the 
dynamical calculation to follow we will use the 12 lowest 
SESs to build the relevant MESs. Besides the vacuum 
MES at zero energy we thus recognize the 12 SESs again 
as MESs occupied by only one electron each. Above the 
energy of the highest SES accounted for we have a rela- 



In a recent papei^ among other things we demon- 
strated a dynamical Coulomb blocking effect in a small 
system where all the relevant SESs are extended. There, 
the mutual Coulomb interaction between the electrons in 
the system prevents the entrance of further electrons until 
the bias is high enough and the occupation of the system 
in the steady state regime shows the well known Coulomb 
steps as a function of the bias between the left and right 
leads. The time-dependence of the contact functions x l 
is described by 

x L>R w = (i - ^4^) > 
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with a 1 = 1.0 ps 1 . We fix the temperature of the reser- 
voirs at T = 0.5 K. 



A. Enhanced occupation by the Coulomb 
interaction 

Here, in a system where not all the states of the 
system are extended we will show that the Coulomb 
interaction can bring about a totally different dynam- 
ical effect: Fig. [4] shows that even for a small bias 
A/i = hl — fiR = 0.1 meV only one electron seems to 
be able to enter the initially empty system in the ab- 
sence of the Coulomb interaction, but the Coulomb in- 
teraction seems to facilitate the entrance of the second 
electron into the system. A glance at Figures [5] and 
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FIG. 4: (Color online) The total charge of the non-interacting 
(0) and the interacting system (I) as function of time. A// = 
0.1 meV, B 1.0 T, V g 1.0 meV, L x 300 nm, and 
gbciw 2 = 30 meV. 

[6] showing the occupation of the MESs \fi) in case of 
the interacting and the non-interacting system, respec- 
tively, indicates a very different charging effect for the 
two cases. In the case of the non-interacting system the 
one-electron MES |5) = |000100000000) is occupied up to 
78% after t = 380 ps and the state |9) = |000000001000) 
carries 12% with the rest distributed to several states. 
For the interacting system the one-electron state |5) = 
1 000100000000) is again initially occupied with a slight 
occupation of |9) = |129) = |000000010000), but 
soon they loose to the two-electron MES 1 19) and 
1 25) that take over. The state 1 19) is a two-electron 
state with the main contributions from |000110000000), 
1 100010000000), 1001100000000), |010100000000), and 
1 101000000000). 1 25) is also a two-electron MES 
with the main contributions from 1 100000100000), 
1 100001000000), |010000010000), |000100100000), and 

joooioioooooo). 

The energy of the state 1 19) is 3.36 meV and the mean 
energy for the same range of time seen in Fig. [7] is higher 
for the interacting case. We assume the leads are in equi- 
librium before the coupling to the wire at t = at a 



FIG. 5: (Color online) The occupation of the interacting MES 
state \fi) as a function of time. A/x = 0.1 meV, B = 1.0 T, 
V g = 1.0 meV, L x — 300 nm, and g^c^J 2 — 30 meV. 
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FIG. 6: (Color online) The occupation of the non-interacting 
MES state as a function of time. A/i = 0.1 meV, V g = 1.0 
meV, B — 1.0 T, L x — 300 nm, and g l al /2 30 meV. 

temperature of 0.5 K corresponding to 0.043 meV. The 
energy 3.36 meV of 1 19) is valid for the case of exactly 
two electrons in the state, here we can only explain the 
occupation of this state well above [1l with the fact that 
it is only partially occupied. In other words the SES ele- 
ments which build the MES 1 19) are only partially occu- 
pied. This corresponds well with the values of the mean 
energy in Fig. [7[ 

In the non-interacting case one electron can enter the 
system and it occupies the state |5) = |000100000000) 
just below jiR. This state is the lowest state with high 
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FIG. 7: (Color online) The mean energy of the non-interacting 
and interacting system as function of time. A/i = 0.1 meV, 
B — 1.0 T, V g = 1.0 meV, L x 300 nm, and g l Q aJ 2 30 
meV. 



FIG. 8: (Color online) The current in the left and right leads 
as a function of time for the non-interacting and the inter- 
acting system. A/i = 0.1 meV, B = 1.0 T, V g = 1.0 meV, 
L x — 300 nm, and gbdw 2 — 30 meV. 



weight in the contact region, below which only states 
(a = 1, a = 2) exist with more weight away from the 
region of contacts, in the sides of the quantum ring, as 
Fig. [3] confirms. (The state a = 3, or |001000000000), 
is nearly degenerate with |5) but does not participate 
in the transport to large extent for the non-interacting 
case). The Coulomb interaction couples together these 
two different types of states and facilities thus the occu- 
pation of two-electron states with one of the electrons in 
a low energy SES state with poor coupling to the con- 
tacts. Similar phenomena is observed at B = 0.5 T if the 
chemical potential in the right lead is placed in a corre- 
sponding location with respect to the 12 non-interacting 
SESs used in the calculation by varying V g . 

The current in the left and right leads displayed in 
Fig. [8] shows that neither the interacting nor the non- 
interacting systems have reached a steady state in the 
380 ps shown. Both leads are still supplying charge to 
the system, but in the case of the interacting system the 
time constants are clearly longer for the charging process 
that is enhanced by the Coulomb interaction. 

The many-electron charge distributions compared in 
Fig. [9] for the non- interacting and the interacting system 
at t = 380 ps confirm this observation and reminds us 
that the two-electron state in the right panel has a rela- 
tively low interaction energy due to the reduced overlap 
of states with high probability in the contact region and 
states with high probability at the other sides of the ring. 



In Fig. 10 we compare the truncated RDO p^ v for the 
case of the non-interacting and the interacting system at 
t = 380 ps. On the diagonal in the left panel we see 
again as in Fig. [6] that for the non- interacting case only 
state 1 5) has considerable occupation, while the vacuum 
state 1 1) is loosing its initial high value and state |9) is 
gaining some weight. So, only one-electron states are 
occupied here. In the case for the interacting system the 
right panel of Fig. 10 shows as Fig. [5] a strong emergent 



6 
4 
2 


>. 

-2 
-4 
-6 



-2 2 
x/a w 




FIG. 9: (Color online) The total many-electron charge den- 
sity for the non-interacting (left panel) and interacting sys- 
tem (right panel) at t = 380 ps. A/i = 0.1 meV, B = 1.0 T, 
V g = 1.0 meV, L x — 300 nm, and goal/ 2 — 30 meV. 



occupation of |19), a two electron state. Figure [2] revealed 
a small energy gap between the single-electron states and 
the many-electron states of the interacting system. The 
GME-formalism excludes any correlation between MESs 
with a different number of electrons. A manifestation 
of this can be seen in the right panel of Fig. [10] where 
only vanishing off-diagonal elements can be found in the 
upper left and the lower right rectangles correlating one- 
and two-electron states. (The one-electron states being 
ji = 1,2,- •• ,12). For the non-interacting case in the 
left panel of Fig. [To] this separation is not as clear cut 
since there the regions of one- and two electrons states 



overlap slightly. More importantly, Fig. 10 reveals a non- 
vanishing correlation between all two-electron states that 
gain any occupation in the system. Initially for both 
cases only the vacuum state |1) was occupied, so clearly 
the coupling to the leads brings about correlation of the 
electrons in the system, and in addition, the Coulomb 
interaction strongly influences this correlation. 
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FIG. 10: (Color online) The reduced density matrix \p^\°' 36 
for the non-interacting (left panel) and interacting system 
(right panel) at t = 380 ps. A/x = 0.1 meV, B = 1.0 T, 
= 1.0 meV, L x 300 nm, and g l Q aJ 2 30 meV. The 
power 0.36 is chosen to make the smaller off-diagonal elements 
visible on scale needed for the larger diagonal elements. 



B. Current oscillations 

In this subsection we will compare the total current in 
the leads for two different values of the magnetic field, i.e. 
at 0.5 T and 1.0 T, and observe how the current changes 
as the bias is increased. We find that the current exhibits 
smooth oscillations with a period of several picoseconds 
for the higher magnetic field as the bias is increased. Here 
we should mention right away that we are not describing 
the small oscillations seen in Fig. [8] at a shorter time 
scale, that are caused by an interference of the coupling 
to different subbands of the leads. (In the case of a one- 
dimensional lead like has been used in the lattice version 
of the GME-model these oscillations do not appear). 4 

The lowest 32 levels of the many-electron energy spec- 
tra for the two different values of the magnetic field are 
displayed in Fig. [TTJ where a care has been taken to iden- 
tify the one- and two-electron MESs for both cases with 
different symbols. Here, A/i = 0.5 meV. We note that 
for B = 1.0 T the one- and two-electron MESs are sep- 
arated by a small energy gap, but not for 0.5 T. The 
time-dependent occupation probability of the MESs is 

Further analysis of the state 



demonstrated in Fig. 12 



structure of the system is done in Table [T] for B = 1.0 T, 
and in Table [H] for 0.5 T. At B = 1.0 T mainly the one- 
electron state 1 5) is initially occupied with lesser prob- 
ability for 1 9). Without the Coulomb interaction only 
these two states with the same components will gain any 
significant probability of occupation. Similarly, as in pre- 
vious section where we analyzed the occupation of the 
system at a lower bias we see here that the Coulomb in- 
teraction facilitates the coupling to the lower lying states 
and thus increases the probability of the occupation of a 
two-electron state with one of the electrons in a low en- 
ergy SES. The next two-electron state to gain significant 
occupation probability is 1 25) with a bit higher energy 
and a higher likelihood for one of the electrons to be just 
above jiR. 

Very similar picture is seen in Table [H] for the state- 
structure at B = 0.5 T, with the exception that only 
a single one-electron state gains significant occupation 
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FIG. 11: (Color online) The energy spectra for the lowest 32 
many-electron states (MESs) built from the lowest 12 SESs 
compared for two values of the magnetic field. The solid hori- 
zontal yellow line indicates the chemical potential of the right 
lead. Different symbols are used for one- and two electron 



states, hr = 1.6 meV. L x - 
meV, and at B = 0.5 T V Q 



300 nm. At B 
= 1.2 meV 



1.0 T Va = 1.0 



probability and the two-electron states coming in after 
the most probable one, |19), are lower in energy than that 
one. Also, though not shown here the non-interacting 
system at B = 0.5 T will gain a slight probability for 
the occupation of a two-electron states. This is caused 
by the missing gap between the one- and two-electron 



states that we earlier pointed out in Fig. [TT 

Figure [13] compared the total current in the left and 
right leads for B = 1.0 T in the top panel for a non- 
interacting and an interacting system. Both show smooth 
oscillations after the initial transient period, but they are 
clearer for the interacting system. The center panel of 



Fig. 13 demonstrates that the oscillations in the current 
are not changed to any extent by smoothly decoupling 
the system momentarily from the leads around t = 150 



ps. The bottom panel of Fig. 13 indicates that for the 
case of B = 0.5 T the oscillations are either absent or 
too weak to be discernible. Clearly, the system is not in 
its ground state, the coupling to the leads moves it out 
of equilibrium, and a glance at Tables |l| and |ll| observing 
what kind of states are available to the system awakens 
the question if the coupling has generated collective oscil- 
lations as in the case of the closed ring subject to strong 
external perturbation.^ In the case of the system with- 
out an interaction the single electron states would have 
a restoring force from the potential defining the wire and 
the ring which do not have a totally flat bottom, see 
Fig. [I] But here the coupling to the leads is weak, the 



8 




FIG. 12: (Color online) For the interacting system the occu- 
pation of MESs as a function of time for B = 0.5 T (upper 
panel), and B — 1.0 T (lower panel). L x — 300 nm. At 
B = 1.0 T V g = 1.0 meV, and at B 



fiR = 1.6 meV, A/2 = 0.5 meV, and g a 



: 40 meV. 



momentary switch-off of it does not influence the oscil- 
lations, and an inspection of the density confirms that 
the density only shows minute oscillations that we will 
describe below. 

The dynamic evolution is governed by the GME ([9| 
and in Fig. [14] we display the correlation of the two- 
electron state 1 19) that gains the highest occupation 
probability with time in our system, i.e. the off-diagonal 
elements of the RDO, p\§, v . For both values of the mag- 
netic field we see indeed oscillations with comparable 
period as the smooth oscillations in the current, but in 
the case of the lower magnetic field many elements show 
strong oscillations but not in phase. For the higher mag- 
netic field only one or two elements oscillate and one of 
them is clearly stronger. 

We see thus that oscillations of the electron correla- 
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100.0% 
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E 9 = 1.88 meV 
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1000000.010000) 


100.0% 
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19) 


E 19 3.36 meV 
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1 101000.000000) 


11.8% 




11) 


010100.000000) 


9.1% 




13) 


001100.000000) 


11.3% 




18) 


100010.000000) 


20.1% 




25) 


000110.000000) 


35.5% 


2e 


25) 


E 2 5 3.68 meV 






34) 


|100001.000000) 


20.2% 




41) 


000101.000000) 


4.9% 




66) 


100000.100000) 


29.0% 




73) 


000100.100000) 


12.3% 




131) 


010000.010000) 


17.8% 



TABLE I: The most probable interacting MES |/x) at t = 380 
ps and B = 1.0 T together with their strongest components 
\v) of non-interacting MESs. The Fock space representation 
of the states \v) is shown with a period indicating the location 
of the chemical potential jiR in the right lead. A/2 = 0.5 meV, 
and gbdw 2 — 40 meV. 
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00110.0000000) 33.8% 
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|01001.0000000) 


4.6% 
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00101.0000000) 
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2e 


13) E13 = 3.02 meV 
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1 11000.0000000) 


88.2% 
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10100.0000000) 


4.5% 
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00110.0000000) 


6.4% 


2e 


18) E 18 = 3.33 meV 






7) 


|01100.0000000) 


10.4% 




10) 


1 10010.0000000) 73.8% 




37) 


100100.1000000) 


11.3% 



TABLE II: The most probable interacting MES at t = 380 
ps and B — 0.5 T together with their strongest components 
\v) of non-interacting MESs. The Fock space representation 
of the states \v) is shown with a period indicating the location 
of the chemical potential /jr in the right lead. A/2 = 0.5 meV, 
and g l ciw 2 — 40 meV. 



tions are inherent in the GME formalism irrespective of 
the presence of the Coulomb interaction or not. It is a 
part of the correlations forced on various states of the sys- 
tem by the coupling to the leads. Here, we observe that 
the Coulomb interaction further couples different types 
of states in the system. States with weak coupling to 
the leads residing in regions of the system away from the 
contact area with states with higher presence in the con- 
tact area. Furthermore, the magnetic field simplifies the 
energy spectrum of the system such that oscillations in 
the correlation of a single pair of MESs will be dominant 
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FIG. 13: (Color online) The total current in the left and 
right leads as function of time for the non-interacting and 
the interacting system at B = 1.0 T (top panel). The total 
current in the left and right leads for a system with steady 
coupling compared to a system where the coupling is momen- 
tarily switched-off smoothly around t = 150 ps (center panel). 
The total currents for a system in B — 0.5 T (bottom panel). 



A/jl — 0.5 meV, L x — 300 nm, and g a 



z 3/2 



40 meV. 



and thus visible in the total current. 

As was stated before the oscillations in the density 
caused by the oscillation in the electron correlation are 
minute. We thus display in Fig. [15] the "derivative" den- 
sity or induced density defined by n(r,t) — n(r,t — St), 
where we have taken St = 1.5 ps instead of comparing 
always to the density at a certain fixed point in time, the 
reason being that the electron charge is still increasing 





FIG. 14: (Color online) The strength of the off-diagonal ele- 
ment |pi9,i/| as a function of time for B = 0.5 T (left panel), 
and B = 1.0 T (right panel). L x = 300 nm. At B = 1.0 T 
V g = 1.0 meV, and at B = 0.5 T V g = 1.2 meV. g l a 3 J 2 = 30 
meV. 



in the time interval used and we see the density peaks in 
the ring away from the contact area always growing. 




FIG. 15: (Color online) The induced density n(r, t) — n(r, t — 
St) at time points, 252.2 ps (upper left panel), 264.3 ps (upper 
right panel), 273,4 ps (lower left panel), and 282.5 ps (lower 
right panel). B — 1.0 T, and g l al /2 = 40 meV. 

In addition, we see that electron density or charge is 
shifted between the peaks and the contact area thus in- 
fluencing the effective coupling between the system and 
the leads. 

In order to check the stability of the results with re- 
spect to the exact location of the window of chemical 
potential we have repeated the calculations for hr = 2.0 
meV. The results are very similar with respect to the os- 
cillations observed in the current and the occupation of 
the MESs. The states gaining highest occupation are still 
the same but at jiR = 2.0 meV additional MESs with a 
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bit higher energy show up with low occupation that were 
almost empty for [ir = 1.6 meV. At present we do not 
feel confident to increase the height of the bias window 
further due to the fact that we have only included 12 
SESs in the calculation. 



V. SUMMARY 

We have used a time-dependent transport formal- 
ism built on the generalized master equation (GME) 
where the mutual Coulomb interaction between electrons 
is treated within the "exact numerical diagonalization" 
or "configuration interaction" to analyze the transport 
properties of a quantum two-dimensional ring or a short 
2D quantum wire with an embedded ring. The quan- 
tum wire is defined by a parabolic confinement potential 
in the wire plane perpendicular to the transport direc- 
tion. The ends of the quantum wire are hard walls that 
are made transparent to tunneling by a weak coupling to 
leads. The shape of the finite wire and the definition of 
the ring potential (JTrj) leads to a ring system with a po- 
tential that does not have a totally flat bottom, see Fig. 
[T] For this reason we have a system that has SESs that 
are ring states for higher energy, but the lowest states 
can be slightly localized in different part of the system. 
We show that indeed, in this geometry the Coulomb in- 
teraction between the electrons increases the occupation 
of the system by coupling states with different localiza- 
tion properties. The correlation of the states in the sys- 
tem caused by the coupling to the leads is enhanced by 
the Coulomb interaction leading to a behavior that runs 



counter to the usual Coulomb blockade in a simpler ge- 
ometry. Of course the Coulomb blocking mechanism is 
inherent in the interacting system, but by comparing the 
interacting system with the noninteracting one we dis- 
cover finer details in the action of the interaction, details 
that are usually collected under the title: correlation ef- 
fects. 

In addition we find current oscillations that are caused 
by oscillating correlation properties of the electrons in 
the system. These oscillations become visible in higher 
magnetic field due to the reduction of MESs active in 
the transport by the magnetic field. As the correlations 
are caused by the coupling to the leads the oscillations 
are visible in systems without or with Coulomb interac- 
tion between the electrons, but the Coulomb interaction 
influences them through its enhancing of correlations by 
coupling of electronic states. 
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